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Abstract. The Laplace transform approach of [GY] is a celebrated advance in valuing 
Asian options. Its insights are fundamental from both a mathematical and a finan- 
cial perspective. In this paper, we discuss two observations regarding the financial 
relevance of its results. First, we show that the [GY] Laplace transform is not that 
of an Asian option price, as reported in [GY] and other papers. We nonetheless show 
how the [GY] Laplace transform can be used to obtain the price of an Asian option. 
Second, we find that following [GY] these Laplace transfoms are available only if the 
risk-neutral drift is not less than half the squared volatility. Using complex analytic 
techniques, we lift this restriction, thus extending the financial applicability of the 
Laplace transform approach. 



1. Introduction 

Asian options are widely used financial derivatives whose valuation has intrigued fi- 
nance theorists for over a decade. In general, these options provide nonlinear payoffs on 
the arithmetic average of the price of an underlying asset. Our understanding of Asian 
options has evolved as an intricate interplay between theoretical and computational 
perspectives. As early as 1992, Yor [Y] expressed the Black-Scholes prices of Asian 
options as a certain triple integral. While this remains one cornerstone of our theo- 
retical understanding of Asian options, the Laplace transform approach of [GY] had 
far reaching consequences for the way we see this valuation problem today. It has first 
spurred the developement of an impressive array of increasingly sophisticated numer- 
ical inversion techniques, see [FMW] for a recent example. It has also motivated the 
development of analytic inversion techniques. For example, by analytically inverting 
the Laplace transform, [SA] improves on Yor's result by expressing the Black-Scholes 
price of an Asian option as a single integral of the product of two functions. To compute 
this integral, [SE] has derived series and asymptotic expansions. The other principal 
computational alternatives available for valuing Asian options are continued fractions, 
expansion into eigenfunctions, and expansion into orthogonal polynomials. Based on 
Yor's results in [Y], this last alternative has been pursued in [Du]. This work has made 
it possible to establish surprising theoretical links between Asian options and modular 
forms, see [Y3], [SL]. Unfortunately, on the practical side [SL] found computing with 
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these Laguerre series to be prone to serious structural difficulties. Thus, at present, all 
workable non-PDE methods for computing exact Black-Scholes prices of Asian options 
derive from the Laplace transform of [GY]. 

It is widely believed that this Laplace transform is that of an Asian option, as can be 
seen from the original paper [GY] or [FMW] for example. This paper shows that this 
is not the case. Nonetheless, the [GY] Laplace transform can be used to obtain the 
price of an Asian option. A more detailed discussion of these somewhat paradoxical 
findings is deferred to Part L To paraphrase those findings here, the function which 
is Laplace transformed in [GY] is a normalized option value in two variables. Both 
variables depend on time to maturity, so the variables are not independent and in fact 
the second depends on the first. This dependency is not accounted for when in [GY] 
the transform is taken over the first variable. 

Thus, the first achievement of the paper is that we show how the [GY] Laplace transform 
can nonetheless be used to value Asian options. The key idea is to replace the problem 
of valuing the Asian option with the problem of valuing a whole family of auxiliary non- 
Asian options. We can then reconstruct the Black-Scholes price of the Asian option 
using this family as a whole. 

Encouraged by the kind support of Yor, we have extended these results from the zero 
dividend yield situation originally considered in [GY] to the case of constant positive 
dividend yields. Part II gives an exposition of this extension, incorporating Yor's tu- 
torials and suggestions. We count it as the second contribution of this paper to the 

understanding of Asian options. 

In providing this extension, we found that the probabilistic arguments of [GY] are 
sufficient to establish the Laplace transforms only if the risk-neutral drift is not less 
than half the squared volatility. This is due to the still limited knowledge of Bcsscl 
processes with negative indices. Unfortunately, this lower bound on the drift restricts 
the financial applicability of the results. For example, if volatility is 30%, we simply do 
not know the Laplace transforms for the non-Asian options if the difference between 
the riskfree rate and the dividend yield is less than 4.5%. Worse yet, the greater is the 
volatility, the greater is this range of parameters in which we do not have the Laplace 
transforms. Unfortunately, it is precisely due to high volatility that Asian options are 
used in the first place. 

Thus, our third contribution is to lift these restrictions in showing the existence of the 
desired Laplace transforms in Part III. As described in §7, we are thus now able to value 
Asian options for all constellations of the riskfree rate, dividend yield, and volatility. 
To establish our result we change perspective form probabilistic methods to complex 
analysis. We apply analytic continuation which is is a standard method from complex 
analysis. However, the analytic continuation of Laplace transforms in the present setting 
appears to be a highly non-standard way of applying this technique, and it has to be 
handled with caution and care. An overview of the key steps of this extension is in §7, 
while the pertinent background from complex analysis is collected in the appendices. 
The reader who perseveres with this article will hopefully join in our realization that 
valuing the Asian option is akin to valuing non-Asian options on the edge of a knife. 
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Part I 

2. Basic notions 

First we discuss the notions basic for the analysis of Asian options. We work in the 
Black-Scholes framework using the risk-neutral approach to the valuation of contingent 
claims. In this set-up there are two securities. There is a riskless security, a bond, 
whose price grows at the continuously compounding positive interest rate r. There is 
also a risky security, whose price process S is modelled as follows. Consider a complete 
probability space equipped with the standard filtration of a standard Brownian motion 
on the time set [0, oo). On this filtered space, we have the risk neutral measure Q, which 
is a probability measure equivalent to the given one. And then we have a standard 
Q-Brownian motion B such that S is the strong solution of the following stochastic 
differential equation: 

dSt = w ■ Sf dt + a ■ Sf dBt, tG[0,oo). 

The positive constant a is the volatility of S. The specific form of the otherwise arbi- 
trary constant w depends on the nature of the security modelled (eg. stock, currency, 
commodity etc.). For example it is the interest rate if iS is a non-dividend-paying stock. 



Fix any time to and consider the accumulation process J given for any time t by: 

J{t) = Sudu. 

Jto 

The European-style arithmetic-average Asian option written at time to, with maturity 
T, and fixed-strike price K is then the contingent claim on the closed time interval 
from to to T paying 



T-to J [ T-to 

at time T. Recall that points in time are taken to be non-negative real numbers. 



3. Normalized valuation of the Asian option 

In the setting of the previous section, the price Ct of the Asian option at any time t 
between to and T is given as the following risk neutral expectation 



JjT) 
T-to 



-K 



which is conditional on the information ^ available at time t. However, following [GY, 
§3.2], do not focus on this price. As described there in great detail, we instead normalize 
the valuation problem as follows. Multiply and divide the strike by the length of the 
monitoring period, then factor out the reciprocal of this length. Split the integral J{T) 
into two integrals, one of which is deterministic by time t and the other of which is 
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random. Couple the deterministic integral with the new strike. For the random integral, 
restart the Brownian motion driving the underlying at time t, and then using the scaling 
property of Brownian motion, change time to normalize the volatility in the new time 
scale to twohundred percent. The precise result is the factorization: 

I — to 

which reduces the general valuation problem to computing 

the normalized time-t price of the Asian option. To explain the notation, A^'^^ is Yor's 
twohundred percent volatility accumulation process 

Jo 

and the normalized parameters are as follows: 

u=^-l, h=^{T-t), q = kh+q*, 
4 

where ^ ^ 

A; = ^, q*= q*{t) = ^ " (t-to) - ^ S^du^ . 

To interpret these quantities, v is the normalized adjusted interest rate, h is the nor- 
malized time to maturity, which is non-negative, and q is the normalized strike price. 
Notice how q becomes affine linear in the time variable h with coefficients k and q* 
depending only on quantities known at time-f. 



4. The pricing dichotomy 

There is now a dichotomy in computing the normalized time-t price 

of the Asian option according to the normalized strike price q being positive or not. 
Indeed, if q is non-positive, Asian options loose their option feature. Computing their 
normalized time-t prices C^*^) is straightforward: 



A 



(^) 



On applying Fubini's theorem, this last expectation is computed as follows: 



eQ 



A 



2{u+l) ' 

for V any real number. The right hand side is analytic in u with its value dX v = —1 

being equal to h. 

In this way, computing C'^'^^ is reduced to the case where q is positive. As a first step 
in this direction, the aim of the present paper is to develop a certain family of Laplace 
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transforms. The valuation problem then reduces to one of inversion. This inversion has 
been effected analytically in [SA] with an integral for the normalized price as the main 
result. Series and asymptotic expansions for computing this integral are given in [SE]. 



5. Intuition about the Laplace transform 

The first key idea for computing the normalized time-t price 



CM = 



of the Asian option is, roughly speaking, to make the maturity date a variable. This is 
the starting point for the Laplace transform approach of [GY] . Since the mathematical 
background is summarized in §17, we try to provide financial intuition here. 

Making the maturity a variable translates into considering C^^) via h{y) = (a"^/4) • {y—t) 
as a function of the maturity y ranging from t to oo. Since option values increase with 
maturity, it seems reasonable to give lower weight to more distant maturities. Translate 
this intuition by continuously discounting y ^ C^'^\h(y)) with a constant discount rate 
z. There is of course no natural choice for z, so we consider it as a variable. In this way, 
we associate to the function of maturity y i— s- C^^){h{y)) a function of this discount rate 
z. Written in full length, this Laplace transform of C^'') is thus formally given by 



/•OO |- 

i?'(CW)(2)= / e'^^'E^ {Ar-{kx+q*)) 
Jo L 



dx 



on changing variables x = h{y) = • (y—t). 



It is a characteristic feature of the "normalized strike price" of the Asian option that it 
depends on the transformation variable x via q{x) = kx+q*. In [GY], however, options 
are considered whose "normalized strike price" is independent of x, and, as discussed 
in §9, this is crucial for the validity of their Laplace transform computations. So it 
appears initially that their results cannot be used to value Asian options. Luckily, they 
can nevertheless be used for exactly this purpose, and this is explained next. 



6. Valuing Asian options using families of non-Asian options 

The basic idea is as follows. Try to reconstruct the normalized time-i price 

CM = £;^[(4"^-(M+g*))+" 

of the Asian option from a family of auxiliary functions whose single members are 
unrelated to the problem of valuing the Asian option but amenable to the [GY] ana- 
lysis. The picture for this is to think of the graph of the normalized price being fibred 
up: try to put through any point of this graph a curve that is amenable to the analysis 
of [GY] . Given that the time dependency of the normalized strike price is the obstacle 
to applying these results, there are natural candidates for such functions. Force this 
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strike price to be constant, and consider for any real number a the function fGY,a on 
the positive real line that sends any a; > to 

/Gy,a(x) = ^«[(4^)-a)+]. 

They are the prices of certain non-Asian options. Taken individually, they cannot be 
used to value the original Asian option. However, as a whole they allow one to recover 
the normalized time-t price. In the notation of §3 this is made precise in the following 

Key Reduction: If q = kh+q* is positive, computing the normalized time-t price 

of the Asian option reduces to computing all fGY,a with a > 0. More precisely, C^*^) is 

obtained by choosing the function fGY,kh+q* and evaluating it at h. 

A moment's reflection will convince the reader that this is true by construction. In 
conclusion we would like to stress that the basic idea of our valuation technique is to 
value an option using a whole family of auxiliary options. And of course, this works in 
much more general situations. 



7. The Laplace transform of fGY,a 

The key reduction of the preceding §6 shows the way to apply the [GY] Laplace trans- 
form to value Asian options. Moreover, adopting the notation of §3, consider the family 
of all functions fGY,a with a > that send any a; > to 

fGY,a{x) = EQ[{A(^)-a)+]. 

Recall that its single members are unrelated to valuing the Asian option, but as a whole, 
the family allows reconstruction of the normalized time-i price C^'^^ of the Asian option. 
As first step in this direction, try to compute the Laplace transform of any fGY,a 

/•oo 

FGY,a{z)= / e-^''fGY,a{x)dx = JZ'{fGY,a){z). 
JO 

Here the complex number z is to be taken in a half-plane sufficiently deep within the 
right complex half-plane such that the integrals arc finite. The function so obtained 
is analytic. The precise conditions under which these integrals are finite is part of the 
description of these generalized modified Geman-Yor Laplace transforms in our 

Theorem: Setting v = 1a~'^w—\, the integrals FGY,a are finite for any complex number 
z with Re (z) > max{0, 2(i/ + l)}. Let z be any complex number with positive real part 
zq such that zq > 2(z^+l) if v >0 and zq > A if u < 0. Then we have 

= z{z-2{.+l)) 

where on choosing the principal branch of the logarithm 

_ 1 

D,{a,z) = ^ J\-ix''+'l^^{^)dx. 
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Herein is the modified Bessel function with complex order ji, as discussed in [L, 
Chapter 5]. The integral D,y{a, z) can be expressed using the confluent hypergeometric 
function $ discussed in [L, Chapter 9]. The precise result is the following 

Corollary: Let z he any complex number with positive real part zq such that zq > 
2(i/+l) if u >0 and zq > 4 if u < 0. Then we have 

Dy{a, z) 



FgyA^) 



where 



Dy{a,z) 



^(v+A+^{z) 

T{^x{z)+^) 



■ $ 



z{z-2{u+l)) 
iy+A+iJi{z 



v+2-ij.(z) 

, ^ l\ (2a) 



1 

e2H 



on setting ii{z) = V^z+u^. 

Wc again stress that these Laplace transforms are not those of the Asian option price, 
but rather the Laplace transforms of the prices of auxiliary options. To obtain Asian 
option prices, we have to invert these Laplace transforms as indicated in §17. In full 
mathematical generality, analytic inversion has been achieved in [S A] . Numerical inver- 
sions have also been accomplished. For example [FMW] computed the following seven 
cases as reproduced in [Du, Table 7.1] 



Case 


r 


a 


T 


So 


20^"^ 


1 


2% 


10% 


1 


2.0 


0.056 


2 


18'A 


30% 


1 


2.0 


0.219 


3 


1.25% 


25% 


2 


2.0 


0.172 


4 


5% 


50% 


1 


1.9 


0.194 


5 


5% 


50% 


1 


2.0 


0.247 


6 


5% 


50% 


1 


2.1 


0.307 


7 


5% 


50% 


2 


2.0 


0.352 



Table 1 Prices 20^"^ for K = 2.0 and 
to = t = using numerical Laplace inversion 

However, the results of [GY] have been applied literally in computing these values. So 
it is only our two results that give a rigorous base for these computations and justify 
their principal relevance and validity for valuing the Asian option a posteriori. 

Obtaining our two results proceeds in two steps. In a first probabilistic step, the ar- 
guments of [GY] are adapted to compute the modified Geman-Yor transforms FGY,a- 
A new exposition of their argument is in Part II incorporating the tutorials and kind 
suggestions of Yor. Its key idea is to factorize the geometric Brownian motion of the 
underlying over a Bessel process of index v. Pertinent notions are discussed in §8. This 
makes time stochastic in such a way that Yor's accumulation process now takes the 
double role of both a stochastic clock and a control variable for the Asian option. At 
first sight, this appears to complicate the original valuation problem. However, this 
double role of Yor's accumulation process is especially suited to the Laplace transform. 
Indeed, using the fact that the "strike price" a is independent of time makes it possible 
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to reduce the computation of the Laplace transform to the fohowing problem: obtain 
explicit expressions for the Bessel semigroup of index u and for a certain conditional 
expectation involving first passage times of Yor's accumulation process. 

However, such results are available for both concepts only if the index v is non-negative. 
This non-negativity condition translates into the condition that the risk-neutral drift is 
not less than half the squared volatility. If this condition is violated, we do not have the 
Laplace transform. Unfortunately, this places restrictions on the Finance applicability 
of the result as dicussed in the introduction. Their effect can already be seen from the 
fact that v is non-negative only in the first two numerical inversions reproduced in 
Table 1 above. 

In a second analytic step, we thus lift this non-negativity restriction on u in §10 to 
§15. For this we use analytic continuation which is a technique from complex analysis. 
To get an idea how it works, the Theorem and the Corollary both ask for establishing 
the equality of two functions. The trick is to treat as a variable too and to consider 
these as functions of both the transformation variable z and u. In fact, the trick is to 
forget the transformation variable as much as possible and consider them as functions 
in V alone. From that point of view first part of the argument establishes the equality 
of two functions in v on the non— negative real line. 

The problem is to establish an analogous statement for v varying over the negative 
real line, and this is what we do in Part III. The idea is to prove more than that. 
This is because functions become amenable to complex analytic methods only on open 
subsets of the complex plane. Thus the identity theorem of §16.2 which we want to 
apply asserts that two functions on a connected open subset of the complex plane are 
equal if they are analytic and agree on a convergent sequence only. So it is in fact no 
longer possible to stick to real numbers only. As a subset of the complex plane they are 
closed with an empty interior, and we have to find ways to thicken the real line. At this 
stage, however, nothing is known about the functions if u is outside the non-negative 
real line. In particular, it is not known if they exist at all, and this we have to establish 
together with their analyticity properties. 

A number of principal problems come up here that have to be be handled with care 
and caution. For instance, the domains of definition in z of these Laplace transforms 
change with the continuation variable v. We have been able to cope with that only by 
establishing yet another upper bound on the price of the Asia option. 

The actual extension argument of Part III is by a two-front attack. The evaluation of 

I?jy(a, z) as a generalized Weber integral using the confluent hypergeometric function $ 
is classical if v is non-negative and z is real. The point is that the identity so obtained 
remains true for complex numbers z and v in the following situation. The real part of 
z is positive and u varies over a band of a certain height depending on Re {z) which 
is symmetric with respect to the real axis. On this domain, the functions involved are 
analytic in u as explained in §11. Crucial for this is our somewhat delicate square root 
lemma of §10. And at this stage we have the Corollary given the Theorem is valid. 
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The analyticity properties in v of the risk neutral expectation L[x, v) = E'^[{^A^x 
are at the base of the Theorem. Originally defined only for non-negative indices v we 
show that this function extends to the whole complex plane as an analytic function 
in V. This determines the analyticity properties of the Laplace transform F[u){z) = 
i?'(L(-, i')){z) which are discussed in §14. The key idea is to consider for any e > the 
strip in the plane of height 2e symmetric with respect to the real axis, the e-thickened 
real line. On choosing z with Re {z) bigger than four and e sufficiently small we can 
show analyticity oi u ^ F{y){z) for elements u in this e-thickened real line with real 
part smaller than e. This e-thickened half-line contains the open interval (0,e). This 



A 





















Figure 1 The e-thickened half-line 

makes possible the analytic continuation in v of the Laplace transforms. Indeed, as 
discussed in §15, having the Laplace transform identity of the Theorem for u in this 
open interval is sufficient to extend its validity to the e-thickened half-line and so to 
all negative indices u. Drawing a picture of the domains, valuing the Asian option so 
is akin to valuing non- Asian options on the edge of a knife. 

Part II Probabilistic part of the proof 



8. Preliminaries on Bessel processes 

As a preliminary to establishing the Laplace transforms of the preceding section, this 
section collects a number of pertinent facts from the theory of Bessel processes. They 
are patterned after the example of the Bessel processes of integer dimension 5>2 which 
are defined by taking the Euclidian distance from the origin of a Brownian motion in 
dimension 5. Applying the Ito Lemma, their squares are the continuous strong solutions 
of the stochastic differential equation 

dpt = 26dt + dBt , po = 1 

[RY, XI, §1]. These stochastic differential equations make sense for any real number 5 
and have a unique continuous strong solution also if 5 is smaller than two. The latter 
is by definition the square of the Bessel process R^^^ on [0, oo) with index v = 5/2 — 1 
starting at with the value 1. While Bessel processes of non-negative indices u stay 
positive if started with a positive value at time zero, Bessel processes of negative indices 
u develop some pathologies. As explained in [RY, XI, §1], in this case they hit zero. If 
— 1 < < 0, they are thereupon instantenously reflected and never become negative. 
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For v = —1 they continue zero. If u < —1, they even continue negative. This is why 
the following Lamperti identity may be surprising 

Lemma: For the index v any real number, we have the factorization: 



e 



Bt+ut 



for any t > 0, where A^^^ = e'^^^'^'^'^^^dw is Yor's accumulation process. 

For > a proof is given in [Y92, §2] while the general case is now contained as 
exercise XI (1.28), p. 452 in the third edition of [RY]. We are indebted to Yor for this 
and for kindly supplying us with the following argument. 

The trick is to apply the Ito rule to the square Zt of If = exp(z/i+iJf) to get 

Zt = 2{y+l) I Zyjdw + 2 f ZyjdByj. 
Jo Jo 

Time change the process using the inverse function r(f) = inf{n|/Q" Zu, dw > t} to Yor's 
accumulation process to get 

r-rit) 

Z,(f) =2(z/+l)t + 2 / Z^dB^. 
Jo 

To interpret the stochastic integral in this sum, apply the basic time change formalism 
for stochastic processes as in [0, §8.5] to obtain 

Mt) rt 

/ ZyjdBw= I Zyj^T'{w)dWyj 

Jo Jo 

where Wt is defined as the stochastic integral Wt = /J^*^*'' ^/Z^, dB^ and is a Brownian 
motion. Using the inverse function theorem of calculus, the derivative of r is equal to 
the reciprocal of the derivative with respect to time of Yor's accumulation process at 
time w. Hence it is equal to the reciprocal of Z^,. On substitution we so identify the 
time changed process Z as continuous solution to the stochastic differential equation 
for the square of the Bessel process of index u: 

=2(z>+l)t + 2 / ^/z^dW^. 
Jo 

Using the uniqueness of the solution of these stochastic differential equations, the time- 
changed process Z so is the square of the Bessel process of index i/. Reversing the time 
change, this translates into 

To establish the identity of the Lemma we have to take square roots. This is not a 
problem if v is non-negative since then the Bessel process takes non-negative values 
only. It does pose a problem if v is negative. In this case, however, recall that the Bessel 
process starts at time zero with a positive value. Since it is continuous by hypothesis, it 
will stay positive until it first hits zero at time t* > 0. Since the accumulation process 
starts with zero at time zero there is a latest point in time t** , infinity admitted, such 
that A('^) is smaller than t* at all points in time t smaller than t**. Thus we have the 
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required identity at least for all points in time t smaller than t** . Now Yt is never zero. 
Since the processes on both sides of the identity are continuous in time, t** must be 
infinity, and the proof is complete. 



9. The modified Geman-Yor Laplace transform for u > 

This section is the first step in the proof of the integral representation of §7 Theorem 
for the Laplace transform 

POO 

FgyA^) = / e-'"fGY,a{^) dx = ^ifGY,a)iz), 

Jo 

where in the notation of §3 and §6 we have fGY,a{x) = E^[{A^'^ —a)'^], for any positive 
real numbers a and x. Wc explain how under the restriction u = 2a~'^w—l > the 
probabilistic arguments of [GY] apply mutatis mutandis to give the 

Lemma: The assertions o/§7 Theorem are valid if v = 2a~^w—l > 0. 

We are very indebted to Yor for correspondence and discussions about this result, and 
are very grateful for his kind support. In the sequel we want to indicate the key steps 
of the proof following [GY] but trying to incorporate his tutorials and suggestions. 
Hopefully no pitfalls have remained undetected. 

The basic idea is to make time stochastic using the Lamperti identity 

for all positive real numbers w as it has been discussed in the preceding section. Here 
is the Bessel process of index u with i?('^)(0) = 1. On applying this Lamperti 
identity, A^'^^ has the double role of both control variable and stochastic clock. That 
the "strike price" a is independent of time now becomes essential. It makes possible to 
transcribe the condition on the control variable to be bigger than a as the first passage 
time 

7"i^,a = (first passage time of Ai") to level a) 
for the stochastic clock. This is the key idea for obtaining the representation 



fGY,aH = F^ 



- 1 



for all w > 0. Indeed, fix any positive real number x, and consider the process A^^^ at 
X on the set of all events where the passage time r^^a takes values less than or equal 
to X. Break the integral defing A^'^^x) at T^^a- The first summand then is A^"^^ at time 
Tu,a and so is equal to a. In the second summand you want to restart the Brownian 
motion in the exponent of the integrand at Tj/^a. Thus shift the variable of integration 
accordingly. The second integral then is the product of exp(2 • {B (t^ a) + i^Tu a)) times 
at x — Ti/^a') by abuse of language after having applied Strong Markov. This last 
process is such that it is independent of the information at time Ti,^a- Unravelling the 
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definition of Tj/a, the first factor is the square of the Bessel process at time a. Now 
taking the expectation conditional on the information at Tj/^a, we thus get: 



The Q-expectation of A'^'^\w) is (exp(2(i/+l)tt>) — l)/(2(z/+l)), as shown in §4 or using 
[Y, §4] . On substitution the required expression for fGY,a follows. 

At first sight this appears to complicate the problem. However, you will see, it is just 
what is especially suited to the Laplace transform FGY,a of fGY,a- 

„2{l/+l)[«;-T..a]+ _ -1 

■ {R^:^f 



FGY,a{z) 



dw . 



Still, for computing this integral you want to interchange the Laplace integral with the 
expectation E^. If z is real it seems best to follow Yor's proposal for justifying this. 
Indeed with the integrand of the double integral in question positive and measurable, 
apply Tonelli's theorem now but justify only in a later step that any of the resulting 
integrals is finite. The case of a general argument z is reduced to this case considering 
the absolute value of the integrand, and the result is the identity 

of measurable functions for any complex number z with Re (z) > 2(i/+l). To identify 
the expectation in the numerator as D,^{a,z) write it out as 



Jo 



■ E^ 



■Pu,a{^,x)dx, 



where pi^^a is the Bessel semigroup of index v starting at 1 at time a. Under the hy- 
pothesis > 0, Yor has in [Y80, Theoreme 4.7, p. 80] (see also [GY, Lemma 2.1 and 
Proposition 2.6]) explicitly computed the conditional expectation factor of this inte- 
grand as the following quotient of modified Bessel functions: 



eQ 



J I I, \a / 



Explicit expressions for the Bessel semigroups p^^a are known only for v > —1, see [Y80, 
(4.3), p. 78] or [GY, Proposition 2.2]. Indeed, in this case the density Pi^^(i(l) f*^) of the 
Bessel semigroup with index v and starting point 1 at time a is 

Pi/,a(l, w) = e • /j. — . 

a \aJ 

So these results are available both only if > 0, and this is where the restriction on v 
comes from. On substitution we then get as required 



Duia,z) 



1 

"2a 



— ]ax . 



To complete the Tonelli argument proposed to us by Yor and to complete the proof, 
we have to establish the finiteness of this integral for any fixed complex number z with 
Re{z) > 2{i/ + l). There is a further technical point to be taken care of. Choose the 
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principal branch of the logarithm to define the square root on the complex plane with 
the non-positive real line deleted, whence 



has a positive real part. Integrability now follows using the asymptotic behaviour of the 
Bessel function near the origin and towards infinity in an essential way. Indeed, from 
[L, §5.7] recall that is a continuous function on the positive real line in particular 
whose asymptotic behaviour near zero is 

Thus, with v + S and Re (//) positive, no integrability problems arise for x near the 
origin. Prom [L, §5.11] we have for large arguments 

/.(O-^ forC-oo, 

whence the factor exp(— x^/(2a)) dominates the asymptotic behaviour of the integrand 
of Di, with X to infinity. Thus the integral is finite, completing the computation of the 
modified Geman-Yor Laplace transform of the Lemma. 



Part III Complex analytic part of the proof 



10. A Lemma about complex square roots 

Lifting the non-negativity restriction on i/ of §9 Lemma to the validity of §7 Theorem is 
by analytic continuation and is achieved in §15. As a prerequisite this section establishes 
the following result about the complex square root associated to the principal branch 
of the logarithm 

Lemma: For any complex number z with Re (z) > 2, we have Re {V^z + v"^ ) > |Re {y)\ 
for any complex number v with |Im(^')| < 1. 

For the proof first write out the complex numbers involved by explicitly introducing 
their real and imaginary parts. So let z = ZQ+izi for real numbers Z{i>2 and zi, and 
v = i/Q+ivi with real numbers {vil < 1 and vq. Then we have 

2zW=^o+i2^i=^ 
setting ^0 = '^zq+Vq — vI and = zi+uqUi. Observe that > S+Uq. 

Moreover recall from §16.1 the working of the complex square root. Write ^ as the 
product of the unique point on the unit circle with angle |^| < vr times the stretching 
factor \^\. Then = |^|i/2 exp(ie/2) whence Re(^V2) = {^{^^ cos{9/2). Moreover 
notice that with positive the absolute value of 6 is smaller than 90 degrees. 
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Distinguish if the absolute value of the argument of ^ is smaller than 60 degress or 
not, i.e., if ^ is inside the wedge in the right half-plane between the line y = 2x and 
y = —2x or not. 

A 




Figure 2 The 60°-wcdgc 

First let ^ be not inside the shaded 60°-wcdgc of Figure 2, i.e., let !^ That the 
real part of ^^/^ is bigger than the absolute value of real part of u then follows by a 
series of estimates. With the absolute value of 0/2 smaller than 45 degrees, cos(0/2) is 
bigger than cos(7r/4) = 2~^/^. For the real part of the square root of ^ we thus get 

Re {V^) = VW- cos Q > ^{il+m" . 

Using ^l<(,l, we have Co+4^i > ^^h whence 

V2^ 



Recalling > 3+fQ, we have > |fo|. Since moreover 5^/^ is bigger than 2-*^/^, 

\Re{v)\, 



as was to be shown. 

Now let ^ be inside the shaded 60°-wedge of Figure 2, i.e., let < ^o- Again we have 
to show that the real part of ^^/^ is bigger than the absolute value of real part of v. 
For this develop the square root at ^ into its Taylor series around 



2i§i_ 



+ 



2' A 3/2 

^- 4 • r]a 



mi? 



with rja = ^o+Q;2i^i for an a in [0, 1]. Then ija is also inside the 60°-wedge. Thus the 
absolute value of the angle 6a of rja is less than 60 degrees. Hence the abolute value of 
39a/2 is less than 90 degrees. Thus cos{S9a/2) is positive, and we have 

1 \ Re(e-*3«"/2-) 



Re 



3/2 
Vol 



Using this non-negativity result in the above Taylor expansion, 

1 



> 0. 



3/2 
Vol 



> V6- 
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Using that is bigger than the square root of S+Vq and thus bigger than |fo|, this 
completes the proof of the Lemma. 



11. First step of the analytic continuation 

As a first step in the analytic continuation argument this section studies the analytic 
properties of the generalized first Weber integral D,y{a,z) of §7 Theorem. On choosing 
the square root associated to the principal branch of the logarithm recall 

for any positive real number a and for any complex numbers z and v such that the 
real part of i^+4+(22:+z^^)^/^ is positive. Taking z with Re > 2(z^+l), we know from 
§9 Lemma that this integral is finite if u is non-negative. The key point is that it stays 
finite if v varies in the band of height two which is symmetric with respect to the real 
axis. Using the confluent hypergeometric function $ discussed in [L, §9.9], the precise 
result is the following 

Lemma: For any complex number z with Re (z) > 2 and any positive real number a, 
we have 

/or any complex number v such that |Im(i/)| < 1 and with ji = {2z+v'^Yf'^ . 

Corollary: For any complex number z with Re (2;) > 2 and any positive real number 
a, sending v to Di,{a,z) defines an analytic map on the set of all complex numbers u 
with |Im {y)\ < 1. 

Proof of the Lemma: For proving the Lemma first establish finiteness of D^{a,z). 
Using the asymptotic behaviour of the modified Bessel function factor of its integrand 
as in §9, it will be finite if the real part of iJ,{z) = (22:+z/^)^/^ is bigger than the absolute 
value of the real part of This is implied by §10 Lemma if the imaginary part of v 

is less than or equal to one and the real part of z is greater than or equal to two. 

For computing this integral we modify the discussion in [W, §13.3, pp.393f] of Han- 
kel's generalization of Weber's first integral. The idea is to expand the modified Bessel 
function in the integrand of 

into its scries of [L, §5.7] and integrate term by term. Using [L, §9.9] this is justified by 
the absolute convergence of the series for the confluent hypergeometric series which is 
to result, and we get 

(2a)M^r(/x+l+n) n! 
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Changing variables y = (2a) ^x^, compute any n-th integral as 

Sorting out the series of the pertinent confluent hyper geometric function we thus get 

Multiplying this expression with exp(— (2a)^^)/a, the Lemma follows. 

Proof of the Corollary: For proving the Corollary, using §10 Lemma notice that /i(z^) = 
{2z+i'^)^^'^ is analytic at any u with |Im(i/)| < 1. Thus the confluent hyper geometric 
function of the Lemma divided by r(/x(^') + 1) is analytic at any such u as composition 
of analytic functions [L, §9.9]. Again using §10 Lemma, the real part of + is 
moreover positive for any i/ under consideration. With the gamma function analytic on 
the right half-plane, thus r((z/ + /x(z/) + 4)/2) is analytic in u, too, as the composition 
of analytic functions. The proof of the Corollary is complete. 

12. Analyticity properties of the moments of A^'") 

This section considers the first two moments £'*^[(A^'' ■')"'] of Yor's accumulation pro- 
cesses Ai") at any positive real number x as functions en,x in the complex variable u, 
and provides the prototype for the discussion to follow. While explicit expressions for 
these moments are in [Y, §4], representations better suited for analytic continuation 
are obtained by direct re-computation. For any first moment function ei^x we have 

^e2x{u+i)_i^ ^ (2a;(z^+l))'^ 
" 2(i/+l) ~''^n=o (n+1)! ' 

originally only if v is different from —1. However, the series is absolutely convergent for 
all v and converges uniformly for i/ in any compact subset of the complex plane. Thus 
ei^x is analytic in i/ on the whole complex plane with ei_a;(— 1) = x. 

Any second moment function e2,x is given by 

^'""^ ' 2{v+l) 2(i/+3) 2{v+l) A{v+2) 

for any complex number v different from —1, —2, and —3. It is an analytic function on 
the whole complex plane whose extensions over the bad points are given at = — 1 by 

g4x(z.+2) _ 
69 r(z^) = 2 ; , 

'""^ ' A{y+2) ' 

at u = —2 by 

e2,x{'') = o/,. I 



2{u+l) 2(z/+3) 2(i/+l) ' 
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and finally at = — 3 by 



9^ , , 9 „4a;(z/+2) 



2{v^\) 2(i/+l) 4(j/+2) 

13. Existence and analyticity properties of the functions jcY.a 

In this section any of the auxiliary functions fcv^a of §6 is considered as function in the 
variable v 

for any fixed positive real numbers a and x. Using §9 Lemma, we know it is defined for 
non-negative real numbers v. However, this has been achieved in a very indirect way 
only: for these values of v the Laplace transforms of the corresponding functions foY.a 
have been shown to be finite. Now more is true indeed 

Lemma: The function L extends to a function on the complex plane which is analytic 
at each point, and for which we have the majorization 



For the proof of the Lemma now set f{w) = {w — a)'^. Applying Girsanov such that 
Wx = ux + Bx becomes a standard Brownian motion, and, dropping reference to this 
new measure, we get 



e ^'^ . 



For establishing the analyticity statement of the Lemma it is thus sufficient to show that 
the expectation factor is analytic in any complex number v. This is true by definition 
if we have the convergent series 



E 



^-^^ ml 

m=0 



f{AP)W^' 



for all ly. For this it is sufficient to show that the scries is absolutely convergent for all 
v. Using the Cauchy-Schwarz inequality this is implied by the convergence of 



m=0 



E 



1U2™ 



for all V. Herein E[f'^{A^x^)] is majorized by the second moment of Yor's zero drift 
accumulation process at x, and so is finite. Since E[W^"^] is majorized by tt"^/^ • 
{2x)"^ ■ m\ for all m > 0, convergence follows using the ratio test. Actually, we have so 
established yet another upper bound to the price of the Asian option. 

To establish the majorization of the Lemma, taking absolute values inside the expec- 
tation in the above Girsanov representation of L gives: 



\L{u)\<E\f{AP).\e 



,2, 



18 



P. Carr and M. Schroder 



The absolute value of the exponential factors are the exponentials of the real parts 
of the respective arguments. Majorizing the function of Yor's zero drift accumulation 
process by this process itself, we get 



Reversing the Girsanov then completes the proof of the Lemma. 



14. Existence and analyticity properties of the functions FGY,a 

While the preceding section has studied the expectations 

L{x,i^) = E'^[{Ai^^-af] 

for any fixed positive real numbers a and x as function in the complex variable v only, 
this section moreover treats a; as a variable. If u is any non-negative real number, we 
have from §9 Lemma that the integrals of the Laplace transforms 

l-OO 

F{v){z) = / e-^^L(.x, ;/) (i.x 
.70 

are finite if Re {z) > 2{u+l). In this section we give independent proofs of the following 
two more general statements 

Lemma: For any complex number v, the Laplace transform F{i')[z) is finite for any 
complex number z with Re (z) > max{0,Im^(z^)/2 + 2(Re (i^) + l)}. 

Corollary: For any complex number z with a positive real part, the map sending v to 
F{i'){z) is analytic in all complex numbers v with Re(2;) > Im^(jv)/2+2(Re {v) + l). 

Proof of the Lemma: As first step in proving the Lemma, we establish for any complex 
number v the majorization 

Jo 

for any complex number z in the sense of measurable functions. Indeed, majorize the 
absolute value of F(i'){z) by taking the absolute value inside the defining integral. The 
absolute value of the exponential factor then is equal to exp(— Re {z)x). Majorizing the 
absolute value of L{x,u) using §13 Lemma, the estimate follows. 

Setting uq = Re (v), now let Rc (z) be positive and such that = Re (z) — Im (j^)^/2 is 
bigger than 2(i/o + l). We compute the Laplace transform of the right hand side of the 
above inequality using §12. If vq is different from minus one, 

poo -I 

using that is bigger than 2{uq+1) to compute the improper integrals. It converges to 
^^q'^ with i^Q going to minus one. Thus it is seen to coincide with the Laplace transform 
for the case vq = —1. The Laplace transforms are finite if is positive and bigger than 
2(fo+l), and then F{i'){z) is finite a forteriori. The proof of the Lemma is complete. 
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Proof of the Corollary: There are a number of reductions for the proof of the Corollary. 
First, it is sufficient to show that F{u){z) is one-fold partially complex differentiable 

with respect to u. For this we can assume that u belongs to any sufficiently small 
bounded open subset U of the complex plane. We claim that computing this partial 
derivative then is by differentiation under the integral defining F{u){z). This will com- 
plete the proof since the integrand of F{i'){z) is analytic in using §13 Lemma. Thus 
we are finally reduced to show that the restriction of the partial derivative of the inte- 
grand of F{iy)(z) to U times the positive real line can be majorized in absolute value 
by an integrable function on the positive real line. 

As a first step, we establish for the expectation L{x, v) the majorization 

<\v\x- eQ [4^^ (-))] + ^/eQ[b^ ■ ^/eQ[{a^^^^ . 

for any positive real x and any v in U. 



{x,u] 



Indeed, apply Girsanov such that Wx = ux+B^ becomes a standard Brownian motion. 
Dropping reference to this new measure. 



L{x,v) = E\f{Af))e 



e 



with f{w) = (w — a)+. The power series in u of the expectation factor derived in §13 

is uniformly convergent on compact subsets of the complex plane. On U differentiation 
of this series with respect to v so is term by term yielding 



E 



The partial derivative of L with respect to v thus equals 



du 



{x,v)=E\f{Af'>)Wxe 



2: ,.2 



XV ■ L{x, v) 



on reversing the Girsanov in the second summand. To estimate the absolute value of 
the second summand use the majorization of §13 Lemma. To estimate the absolute 
value of the expectation in the first summand, apply Cauchy-Schwarz to get 



E 



/(4°))W^,e'^^--f'^1 < E[W^] -eI /(40))e''^'"-f'^' 



2i 



Majorize the second expectation on the right hand side arguing as for §13 Lemma, and 
the desired inequality follows. 



As a next step, we show that after shrinking U if necessary the absolute values of the 
functions x i-^ ex.p(—zx){dL/du)(x,u) with u inU can be majorized in dependency on 
[/ by a function which is integrable on the positive real line. 
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For this apply the inequahty of the first step. There we have E'^ [i?^] = x, for any 
positive real number x. Using §12 for the other expectations, for any u in the 
absolute value of the partial derivatives x ^ {dL/du){x,u) can thus be majorized by 
a function on the positive real line with the following properties. It is continuous 
in both the parameter Re (n) and its argument x, it behaves like a polynomial for 
small values of x, and for large values of x behaves like exp(mux) on abbreviating 
ruu = max{0, Im ^(u)/2+2(Re (n)+2)}. Thus x i-^ exp(— zx) • {dL/du){x, u) is intcgrable 
on the positive real line if Re (z) is positive and bigger than Im^(7i)/2 + 2(Rc {u) + 2). 

Now recall we had fixed any v such that Re (z) > Im^(z/)/2 + 2(Re (z/) + 2) and then 
chosen a bounded open neighborhood U of v. Shrink U if necessary so that Rc {z) > 
Im^(M)/2 + 2(Re (w) + 2) for any u in its closure. Choose u* in this closure with rriu* 
maximal. This takes care of the growth behaviour of the exponential factors in the 
fu in such a way that the ones in f^* dominate those of all other functions for 
larger values of x. Now take stock of the functions that occur in /„* as coefficients of 
the exponentials. They are continuous in Re (n) in particular. Letting u vary over the 
compact closure of U majorize them by constants also depending on U so that the 
resulting function fu majorizes all other functions fu with u in U. Since fu behaves 
for large values of the argument x like cxp(mu.x) too, x i— >• exp(— xRe (z)) • fu{x) is 
integrable on the postive real line by the choice of u* . So it qualifies as a majorizing 
function we have been looking for, and this completes the proof of the Corollary. 



15. A summing up 

It remains to pull things together and establish the two results of §7. 

First notice that §7 Corollary is implied by §7 Theorem using the computation of 
D^{a, z) in §11 Lemma. Thus we are reduced to prove §7 Theorem. 

The proof of §7 Theorem is by analytic continuation using the identity theorem recalled 
in §16.2. As a consequence of §9 Lemma it remains to establish this result for negative 
indices u only. The existence of the Laplace transform on all complex numbers z with 
positive real part required in the Theorem then has been proved in §14 Lemma. 

For establishing the crucial Laplace transform identity of §7 Theorem, thus let z be any 
complex number with real part Re(z) > 4 and choose < e < 1. Using §11 Corollary, 
the generalized first Weber integral 

DAa,z) = '^J\-i,^^^I^(^dx 

is analytic in u on the e-thickened real line which consists of all complex numbers 

v with |Im < £. Picture as the band of height 2e symmetric with respect to the 
real axis. If we choose e so small that Re {z) > 2(2+2£), we claim to have analyticity 
of the Laplace transform 

/•oo 

F{i^){z) = / e-^''L{x,u)dx 
Jo 
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as a function in v on the subset of which consists of all complex numbers v with 
|Im {v)\<e and Re {v) < e. Indeed, if Re {u) < e, we have 2(2+2£) > 2£+2(Re (i^)+2). If 
|Im (z^)| < e, wc have 2e > Im^(i/)/2 since e < 1. Any u in so satisfies the inequality 
Re (z) > Im^(i/)/2+2(Re (z^)+2) of §14 Corollary, and the claim follows. 

So the functions F(i'){z) and 

^>'"^ = z(.-2(.+ l)) 

are analytic as functions in u on the £-thickened half-line B^. It is now a consequence 
of §10 Lemma that we have F{ij){z) = D*{a, z) for all > such that 2(i/+l) < Re {z). 
With Re (z) > 4 this so holds a forteriori for all > such that 2(z^ + l) < 4, i.e., 
for all non-negative real numbers v smaller than 1. With e smaller than 1 we so have 
F{v){z) = D*{a, z) for any v in the subinterval (0, e) of B^. With B^ open and connected 
the identity theorem literally applies to give that F{v){z) = D*{a,z) for all i' in Sg. 
This identity then holds a forteriori for all real numbers v in B^, i.e., for all v < e. The 
proof of §7 Theorem is complete. 



16. Appendix: notions from complex analysis 

This section reviews pertinent results about complex analysis. It is not meant as a 
substitute for a closer study of basic texts like [FB] , [C] , or [R] . 

On choosing any root z of a;^+l, any complex number z has the form z = x+iy where 
the real number x = Re (z) is its real and the real number y = Im (z) its imaginary 
part. The complex conjugate of z has the same real part as z, its imaginary part is the 
negative of that of z. The product of z and its complex conjugate is a non-negative 
real number whose square root is the absolute value \z\ of z. Thinking of z as the pair 
with entries the real and imaginary part, it equals the Euclidian distance of z to the 
origin in the plane. 

16.1. Exponential function, logarithms, and square roots 

In analogy to analysis over the reals we have the notion of a complex-valued function 
/ being differentiable in an element of some open subset of the complex plane C where 
it is defined. In contrast to real analysis, however, complex differentiability is equiv- 
alent to analyticity in the sense that if / is differentiable in a point of its domain of 
definition, it is already differentiable to any order there and the resulting Taylor series 
converges in a neighborhood of that point. A basic example of an analytic function is 
the exponential function exp given by the series exp(2;) = X^^q ^"'/'^' for any complex 
number z. It extends the exponential function on the real line. To get accustomed to 
its properties, the absolute value of the exponential function at any complex number z 
equals exp(Re {z)), the value of the exponential function at the real part of z: 

Indeed, the square of the absolute value of exp(2:) is the product of exp(2:) with its 
complex conjugate. Using the series, the complex conjugate of exp(2:) is the exponential 
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function at the complex conjugate of z, and the above relation follows. The fundamental 
problem now is that this function is periodic and so has no inverse. Indeed, performing 
the series calculations, 



for any angle 0. There is a way out of this. Think of exp(i^) as the element on the unit 
circle corresponding to the angle 6. Any complex number z different from zero so is a 
product of its distance \z\ from the origin times an element exp(i^), just draw a ray 
through it and the origin to see this. So z = exp{i9 + log \z\) which forces 



The problem is that 9 is determined only up to integer multiples of 27r. To get rid of this 
indeterminancy put restrictions on the range of values of the angles. In the sequel the 
convention 16*1 < tt is chosen. The upshot is that so you get a continuous logarithm not 
on the whole complex plane but only on C_ = C \ R<o, the complex plane cut along 
non-positive real line. This is the principal branch of the logarithm. For any element z 
there is now a unique angle 9 = arg(2:) of absolute value less than tt, the argument of 
such that z = \z\ exp(i0). This yields the square root: 



The effect of this square root associated to the principal branch of the logarithm can 
be pictured as follows. Any element in C_ is on a unique ray through the origin with 
angle 9 between minus tt and tt in a distance R from the origin. The effect of taking its 
square root is then halving the angle of this ray, so putting things into the right half 
plane, and changing on this new ray the distance of the image under the square root 
to the origin to the square root of the old distance R [FB, I§2]. 

In particular notice what happens at any point of the cut. The limit from the upper 
edge of the cut to this point of the square root differs from that taken from the lower 
edge of the cut by a factor minus one. Indeed, the angles corresponding to the upper 
and the lower edge of the cut differ by 27r which upon taking square roots produces the 
factor exp(z7r) and thus gives the minus one by which the limits differ. 

16.2. The identity theorem 

The identity theorem asserts that any two functions analytic on an open and connected 
subset G of C are identical if they agree on a sequence in G which converges in G. 
Indeed, it is sufficient to prove that any function / as in the theorem is identical to 
zero if it vanishes on any convergent sequence dn as in the theorem. Suppressing some 
purely topological nicities, the connectedness of G implies that it is sufficient to prove 
/ identical zero in a on a neighborhood of the limit d of this sequence. So develop / 
in its Taylor series around d. This Taylor series converges on an open neighborhood B 
of d in G. Since / is continuous and vanishes on the points dn, it is zero at d whence 
the zeroth above Taylor coefficient is zero. Dividing f hj z — d so gives a function on 
B which is analytic there. It equals zero on all d„ different from d, and so is zero on 
d, too. Consequently the first Taylor coefficient of / is zero. Proceeding inductively, all 
coefficients in the above Taylor expansion of / are zero. Thus / is zero on B. 



exjp{i9) = cos{9)+ism{9), 



log(z) =i9 + log \ z 




arg(z) 
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17. Appendix: Laplace transform 

This section collects a number of relevant facts about the Laplace transform from [B], 
[D] . The class of functions considered in the sequel is that of functions of exponential 
type, i.e., of continuous, real-valued functions / on the non-negative real line such that 
there is a real number a for which exp(at)/(t) is bounded for any t > 0. The Laplace 
transform is the linear operator that associates to any function / of exponential 
type the complex-valued function ^{f), analytic on a suitable complex half-plane, 
which for any complex number z with sufficiently big real part is explicitly given by: 



To fix ideas, consider the function given by f(t) = exp(ai) with a any positive real 
number so that formally 3'{f){z) = e^p{—{z — a)t) dt. Recall from §16.1 that the 
absolute value of cxp(?i;) is cxp(Rc (id)), for any complex number w. This improper 
integral thus exists iff the real part of z — a is positive, i.e., iff the real part of z is 
bigger than a. In this case, S'{f){z) equals (z — a)^^ and is analytic on the half plane 
{2;|Re(2;) > a}. Now mind the following fallacy! The function z i— > {z—a)~^ is defined 
for any complex number z different from a. However, this function has no relations 
whatsoever with the Laplace transform on this half-plane. In fact, with the defining 
integral of the Laplace transform being infinity on {z\Rq {z) < a}, the latter does not 
make sense there at all. 

That i?' is an injection is a consequence of ultimately the Weierstrafi approximation 
theorem. Its inverse, the inverse Laplace transform is expressed as a contour 

integral by the complex inversion formula of Riemann. It applies to any function h 
analytic on the half-plane {Re [z] >Zq} with zq any sufficiently big positive real number. 
Then it asserts that for any positive real number t we have: 



if h moreover satisfies a growth condition at infinity such that the above integral exists. 
For any function / of exponential type we so have in particular o ^){f) = f. 
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